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Abstract 


Less than catastrophic meteoroid impacts over 10 5 years may change the shape of small 
rubble-pile satellites in binary NEAs, lengthening the average BYORP (binary 
Yarkovsky-Radzievskii-Paddack) rate of orbital evolution. An estimate of shape-shifting 
meteoroid fluxes give numbers close enough to causing random walks in the semimajor 
axis of binary systems to warrant further investigation. 
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1. Introduction 

Perhaps ~15% of the near-Earth asteroids (NEAs) are binaries (Pravec et al. , 
2006). The dynamical lifetimes of these small binaries is thought to be ~10 5 yr due to the 
binary Yarkovsky-O’Keefe-Radziebskii Paddack (BYORP) effect (Cuk and Burns, 
2005). BYORP either increases the system’s semimajor axis to the point where the 
satellite escapes, or decreases it until the satellite reaches the primary. The large 
percentage of binaries may mean that binaries are created at a high rate, destroyed at a 
slow rate, or the dynamical lifetimes are longer than ~10 5 yr, or some combination of 
these. 

One reason for believing the dynamical lifetimes are longer is that the 10 s years 
refers to uninterrupted BYORP evolution. If occasionally a synchronous satellite were to 
have an episode where the satellite flipped and the other end thereafter synchronously 
pointed at the primary, then the BYORP effect would change sign and the evolution 
would be reversed until the next flip, lengthening the dynamical lifetime via a random 
walk. There purely dynamical reasons why such end-for-end switching might occur (Cuk 
and Nesvorny, 2010). Also, a satellite may accrete mass from the primary; 1999 KW4 
might speed up from the YORP effect, for example (Ostro et al., 1996), and throw off 
pieces which end up on the satellite, changing its shape and altering BYORP. 

We examine here another possible mechanism as to why steady BYORP 
evolution might be interrupted: less than catastrophic meteoroid collisions with the 
satellite of an NEA. The idea is this: most asteroids appear to be rubble-piles (e.g., 
Harris et al., 2009). BYORP depends on the satellite’s shape. If the satellite is a nibble- 
pile, then a collision may disturb the nibble structure enough to change the shape of the 
satellite, altering the BYORP effect, including perhaps changing its sign. If the collisions 
happen often enough, then once again a sort of random walk will ensue, lengthening the 
average dynamical lifetime. The amount of energy needed for an impactor to 
significantly modify the shape of small rubble-pile satellites is estimated below, along 
with the frequency of such hits. Changes in satellite spin state are ignored. Detailed 
modeling of the impact process is beyond the scope of this paper. 
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2. Gravitational potential energy 

We forsake detailed computations of shape changes here for a simple estimate 
based on the energy required to pull apart a spherical rubble-pile held together only by 
gravity; i.e., a self-gravitating sphere composed of millions of frictionless, non-sticky 
marbles. The gravitational potential energy U of a sphere of radius R and uniform density 
p is 
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3GM 2 3G/4jt 
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48 
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( 1 ) 


where G is the universal constant of gravitation (e.g., Rubincam, 1979). This is the 
energy required to pull apart all the marbles composing the sphere and send them to 
infinity, where they arrive with zero velocity. Let 

AU = | U\F (2) 

where AU is the change in energy necessary to significantly alter the BYORP effect, 
such as cancel, double, or reverse the sign of the effect.; thus AU is expressed as a 
fraction F of U in the above equation. 

The value for F is estimated via an energy argument. Consider a non-rotating 
spherical satellite of mass M and radius R as above. Suppose an impact is energetic 
enough to fission the satellite into a contact binary consisting of two non-rotating 
spherical asteroids as shown in Fig. 1. Neither the assumed original spherical satellite nor 
the resulting contact binary will actually exhibit BYORP; they are too symmetrical. The 
purpose of the estimate is to make an educated guess based on easily computed shapes, 
without actually constructing an explicit model. 

In the contact binary, assume the larger object has mass M\ and radius R\, and the 
smaller has mass M 2 and radius Ri = JR. Being a contact binary, Mi simply rests on M\$ 
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surface. Also, both M\ and Mi have the same uniform density p as M. Moreover, assume 
the total mass and volume are conserved. In this case i ? 3 = R\ 2 + Ri, so that R\ = (1 -/ 
3 ) l/3 /?, Mi =f 2 M, and M\ = (1 - f i )M. The gravitational potential energy is now 

v ,_ 3 GM] 3 GM\ GM ] M 2 

5 /?, 5 R 2 + R 2 1 


where the last term is the interaction energy of the two objects. The factor F as a function 
of /is then 


U'-U 
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and is shown in Fig. 2. An educated guess is that / « 0.315, so that F = 0.01; thus the 
gravitational potential energy changes by 1% for this value of f 

3. Impactor mass and velocity 

When an impactor hits the satellite, it is the explosion that provides the energy for 
the shape change. It is necessary to estimate the impactor’ s mass in order to compute the 
flux of impactors, because the flux depends on mass. The impactor’ s kinetic energy is 
mv 2 /2, where m is its mass, and v is its velocity. Let 0 < £ < 1 be the efficiency with 
which the impactor’s kinetic energy is transformed into the required potential energy A U. 
In this case 



-Mr- 


1 

5 R 
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Solving for m yields 
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so that the lower the efficiency £, the greater the impactor mass needed to effect shape 
changes. 

Throughout the following it will be assumed that the satellite has p = 2000 kg m' 3 . 
The low density p is due to the voids in the rubble-pile. 

To get an idea of the magnitude of m in (6), it is necessary to estimate v, F, and §. 
Velocity v is the least uncertain of these. Suggs et al. (2010) find that objects in the 0.1 - 
10 kg range impacting the Moon overwhelmingly come from meteoroid streams in 
comet-like orbits. It will be assumed here that impactors slightly smaller than 0.1 kg also 
follow these orbits. The Taurids have relatively low velocities of -2.4 x 10 4 m s' 1 , while 
the Leonids have extremely high velocities of ~7.1 x 10 4 m s' 1 (e.g., Moser et al., 2011, 
Table 3; Kokhirova and Borovicka, 2011). The Geminids have intermediate velocities 
of -3.5 x 10 4 m s' 1 (e.g., Moser et al., 2011). A Geminid-like value of v = 3.333 x 10 4 m 
s' 1 is adopted here for v (so that conveniently v 2 = 10 9 m 2 s' 2 ). 

Of the list of 37 near-Earth binaries compiled by W. R. Johnston, many have 
satellites with R * 150 m, while the largest on the list, 2001 SN 263 , has R = 500 m 
(http://www.johnstonsarchive.net/astro/astcroidmoons.html). For a satellite with R = 100 
m, M= 8.4 x 10 9 kg, V esc = (2 GM/R)' 12 = 0.1 1 m s' 1 , g = GM/R 2 = 5.6 x 10' 5 m s' 2 , and U 
= -(3/5) GM 2 /R = -2.8 x 1 0 7 J, where V esc is the escape velocity and g is the surface 
acceleration of gravity. For such a 100 m satellite, 

m — 0.056 (F/£) kg (7) 

by (6). For the likely values of F < 0.1 and £ > 0.001, impactor masses fall in or a little 
below the range studied by Suggs et al. (2010); so choosing v as a meteoroid stream 
velocity seems self-consistent and reasonable. 

The factor/ 7 has already been estimated to be -0.01. It remains to guess £. Kinetic 
energy may be lost through the escape of ejecta. An impact on a satellite made up of only 
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regolith would be able to scatter loose material, but the medium is highly dissipative. An 
impact on a rubble-pile composed of boulders would expend much energy in breaking the 
chemical bonds holding the rock together, but the seismic shaking might jostle the 
rubble-pile enough into changing shape (e.g., Sanchez and Scheeres, 2009). If the 
process is only 1% efficient regardless of composition, then £ = 0.01, which gives F/% = 
1 . Plugging this in (7) gives for the 100 meter satellite m = 0.056 kg, which is a mere two 
English ounces (2 oz.). 

From these considerations it seems sensible to assume that F/% may be 
somewhere between the limits of 0. 1 < F/h, < 10. This range is explored below. 

4. Impact flux 

The impact rate on the satellite will be found next. For the Earth a good 
approximation is (e.g., Bland, 2005) 


N f (> m) = B F m~ a (8) 

where Ne{> m) is the number of meteoroids per year which enter the Earth’s atmosphere 
with mass > m, and Be and a are constants. The above equation can be written 

log Ne (> m ) = -a log m + log Be (9) 

where the logarithm is base 10. For meteoroids with masses in the range 1 kg < m < 10 
kg 

log Ne,\ = -0.926 log m + 4.739 , (10) 

(Bland, 2005; Bland and Artemieva; 2006), so that in this case a = 0.926 and log B E = 
4.739. (See also Silber et al. (2009) for similar parameters.) 

Curiously, there appears to be a break in the curve at smaller masses. Bland 

(2005) and Bland and Artemieva (2006) find that 
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log Ne ,2 (> m) = -0.480 log m +4.568 (11) 

for m < 1 kg, based on the fireball data of Halliday et al (1989). According to (11), 
meteoroids are less numerous than would be expected by extending (10) to masses < 1 

kg- 

The lunar impacts tend to support (11) more than they do (10). Suggs et al. (2011) 
find that Ne is -150,000 yr' 1 (their Fig. 6), based on the lunar data for m > -0.1 kg and 
with no account being taken of the greater focusing factor of the Earth relative to the 
Moon (R. M. Suggs, private communication, 2012). Assuming a focusing factor of 1.11 
for the Earth (see below) and ignoring that of the Moon, their figure rises to -165,000 yr' 
', which is 50% higher than (11), but still closer to the Ne = 1 1 1,686 yr' 1 given by (1 1) 
than the Ne = 462,381 yr" 1 given by (10). The lunar data thus give some support for the 
downturn in flux for masses less than 1 kg. However, it should be noted that the lunar 
lower limit of 0. 1 kg is approximate, and the value of Ne based on lunar data may go up 
or down depending on how accurately the lower limit of 0. 1 kg can be established. 

Assume that the binary system is an NEA and an orbit about the Sun which is 
orbit similar to that of the Earth; then the impact rate N,{>m ) on the satellite will be 
approximately 


Ni(>m) = 4jzR 2 N E ,i(>m)/(oopik A E ) , 


( 12 ) 


which can be written 


N,(> m) = 


4 jtB E 
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using (6), where i = 1 if the impactors follow (10) and i = 2 if they follow (11). Also, Ae 
= 5.1 x 10 14 m is the area of the Earth, 4 nR 2 is the area of the satellite, and 
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is Opik’s (1951) gravitational focusing factor due an object’s gravity, which pulls in 
meteoroids that would otherwise miss it. Here v eS c is the Earth’s escape velocity, and v the 
impactor velocity at a great distance from Earth. Substituting v esc - 1.1 x 10 4 m s' 1 and v 
= 3.3 x 10 4 m s' 1 in the above equation yields ao P ik = 1.111; this value is used here. The 
focusing factor for the satellite is taken as unity because small asteroids have so little 
gravitation. 

If the break in the curve is somehow incorrect and (10) holds for m < 1 kg, then 
Be = 54828 yr' 1 and a = 0.926, so that for the satellite 


NA> m) = 17.5x10“ 


-JW0} 263 ' 


F 


0.926 


yr 


(15) 


On the other hand, the data suggest that (11) more nearly holds. In this case Be = 36982 
yr' 1 and a = 0.480, so that the number of hits per year on the satellite would be 


NJ> m) = 3.3x 10“ 


100 


R 


.0 480 


yr 


(16) 


5. Results 

Let 

Hi = \0 5 N\(>m) (17) 

be the number of hits in 10 5 yr when the meteoroid masses follow N\, and 


H 2 = 10 s N 2 (>m) 


(18) 
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be the number of hits over the same time interval if the meteoroid masses follow N 2 . Hi 
and H\ are plotted in Fig.s 3 and 4, respectively, as functions of R, for the values F/£ = 
0.1, 1, and 10. 

The mass m required for an impactor to substantially alter the shape of a 100 m 
rubble-pile found above is surprisingly small. For the assumed value of v and for F/£ = 1, 
it is only 0.056 kg. However, F/£ = 1 is little more than an educated guess. Can the mass 
really be as small as the amount of loose change in one’s pocket? 

Assume R = 100 m. If the downturn in flux suggested by the terrestrial fireballs is 
indeed real, then (16) holds, and H 2 = ~3 by (18) for F/§ = 1 (see Fig. 3). The shape 
changes would affect BYORP evolution only a little; there would not enough hits over 
10 5 years to give much of a random walk in the semimajor axis of the orbit. However, the 
high (and perhaps unlikely) value F/£ = 10 would lengthen the dynamical lifetime by 
maybe something like a factor of ( H 2 ) u 2 = 10 1 ' = ~3.3, which would be more significant. 

On the other hand, if the dearth of small impactors is somehow bogus and (17) 
gives the real, much higher flux, then the number of hits could be much higher: H\ = 17.5 
for R = 100 m and F/£ = 1 (see Fig. 4). This would lengthen the average dynamical 
lifetime by perhaps a factor of ~(Fi) l/2 = (17.5) 1/2 = ~4 over uninterrupted BYORP, 
which would be substantial. The number of hits II\ can be very high for F/§> 1, giving a 
long random walk (Fig. 4). 

The lunar impact data are in between H\ and H 2 . For instance, the flux is 50% 
greater than that suggested by the fireballs. For R= 100 m and F/£j= 1, the number of hits 
in 10 5 yr is 1.5 x H 2 = 5, still rather small for a random walk. 

It should be noted from (15) that the flux of shape-changing impactors can depend 
strongly on satellite radius R if the meteoroids follow N\; the flux varies like l/i? 2 ' 63 and 
the size of the satellite could not be much larger than R = 300 m to have a significant 
number of hits in 10 5 years, even if F/£= 10 (Fig. 4). On the other hand, if the meteoroids 
follow N i, then (16) holds, and the flux drops off more gently, like 1/P 0 4 ; but the flux is 
already small for R = 100 m, regardless of F/£, and smaller still for satellites with R > 
300 m (Fig. 3). 
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6. Discussion 

The upshot of these results is that the number of shape-changing hits needed to 
produce a long random walk depends sensitively on whether the meteoroid flux follows 
N\ or N 2 ,, and what the value of Fl% is. Most of the satellites listed at W. R. Johnston’s 
website have R < 300 m, so that for them random walks may be possible. However, 
random walks for satellites with R > -300 m are unlikely regardless of the flux or the 
value of F/§; the meteoroid masses are too big to give much of a flux for these larger 
sized satellites. 

It may be that a random walk is not necessary to lengthen the dynamical lifetime. 
A non-equilibrium shape for the satellite is required for BYORP to operate. Perhaps only 
one or a few hits are required to jostle a pile into something like an equilibrium shape 
(Tanga et al., 2009). In other words, the arrow points the wrong way in Fig. 1. A 
perfectly equilibrium shape would shut down BYORP completely and the orbit would 
stop evolving via BYORP. 

The whole argument here assumes that the satellites are rubble-piles. However, 
asteroids smaller than -50 m are often solid objects, as deduced from their rapid rotation. 
Rotationally fission asteroids (e.g., Jacobson and Scheeres, 2011) may simply throw off 
a monolith, which becomes the satellite. Hence the results found above may apply mainly 
to small satellites, but not too small so that they are monoliths. 

So, do shape-changing impacts significantly lengthen the dynamical lifetime of 
binary NEAs? Given the uncertainties in the assumed parameters, nothing conclusive can 
be derived from the above analysis. However, the results are tantalizing enough to 
warrant a closer look. Do the meteoroid masses really follow N\ as given by (15), or Ni 
as given by (16), or something in between, as suggested by the lunar data? What are 
reasonable values for F/§? Granular models such as those of Schwartz et al. (2012); 
Tancredi et al. (2012), Sanchez and Scheeres (2011), and Richardson et al. (2009) 
would be of great aid in settling what might be reasonable values of F/%. 
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Figure captions 
Figure 1 

The spherical satellite (left) fissions into a contact binary (right) after colliding with a 
meteoroid. Neither the original object nor the objects in the contact binary have any 
rotation. The mass and volume of the original satellite are conserved, so that no mass is 
ejected. The dimensions in the figure are approximately correct for the smaller contact 
object having/= 0.315 times the radius of the original satellite. 


Figure 2 

The fractional change F in the gravitational potential energy U after the spherical satellite 
of radius R fissions into a contact binary, plotted as a function of f where the smaller 
contact object has a radius JR. The mass and volume of the original satellite are 
conserved, so that no mass is ejected. The star indicates / = 0.315, the value for which F 
= 0.01. 


Figure 3 

H\, the number of hits in 10 5 years for meteoroids following A/j, as a function of/?. The 
open circles are for ( F/l ;) = 0.1; the solid circles are for (F/f) = 1.0, while the open 
squares are for (F/§) = 1 0. 


Figure 4 

Hi, the number of hits in 10 s years for meteoroids following Ni, as a function of R. The 
open circles are for (F/£) = 0.1; the solid circles are for (F/g) = 1.0, while the open 
squares are for (F/§) =10. 
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